Quantum mechanical relaxation of open quasiperiodic systems 
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We study the time evolution of the survival probability P(t) in open one-dimensional quasiperi- 
odic tight-binding samples of size L, at critical conditions. We show that it decays algebraically as 
P(t) ~ t~ a up to times t* ~ L 7 , where a = 1 — Dq , 7 = 1/Dq and D§ is the fractal dimension 
of the spectrum of the closed system. We verified these results for the Harper model at the metal- 
insulator transition and for Fibonacci lattices. Our predictions should be observable in propagation 
experiments with electrons or classical waves in quasiperiodic superlattices or dielectric multilayers. 
PACS numbers:05.60Gg,03.65Nk,72.15Rn 
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The decay properties of open quantum mechanical sys- 
tems, have been attracting considerable attention for sev- 
eral decades. Their study was motivated by various areas 
of physics, ranging from nuclear jjj, atomic 0] and molec- 
ular H physics, to mesoscopics Q and classical wave 
scattering [6L In recent years, the interest in quantum 
mechanical decay was stirred by mesoscopic cavities and 
microwave billiards where immediate experimental real- 
izations have become feasible 0. At the same time var- 
ious analytical techniques have been developed to study 
the problem in more detail. One possible formulation of 
the problem is to consider the survival probability P(t) 
of a wave packet localized initially inside the open sys- 
tem. In particular it was found that this quantity ex- 
hibits slower than exponential decay (i.e. long-time tails) 
for disordered wires in the localized (in one-dimension) 
and metallic (in higher dimensions) regimes 
Moreover, for ballistic systems, Random Matrix The- 
ory (RMT) predicts an algebraic decay P(t) ~ l/t 0M / 2 , 
where M is the number of channels and (3 — 1(2) for 
preserved (broken) time reversal symmetry 12|l3f|. 

The investigation of the survival probability has re- 
cently been extended to quantum systems with a mixed 
classical phase space [p 15 where it was found that 
P(t) ~ 1/t. The same algebraic decay was found for sys- 
tems with exponential localization. In both cases, this 
law is related to localization and tunneling effects and 
applies for intermediate asymptotic times ]l4| ]. 

The subject of the present paper is the survival prob- 
ability in a new setting, namely a class of systems whose 
closed analogues have fractal spectra. The latter exhibit 
level clustering jl6| in strong contrast to the level repul- 
sion predicted by RMT for systems in the ballistic regime 
and to the Poisson statistics in the localized regime jL7) . 
Representatives of this class are quasi-periodic systems 
with a metal-insulator transition like the Harper model 
pq,|8|, Fibonacci chains ]lq,[l9|, or quantum systems 



with a chaotic classical limit like the kicked Harper model 

Here for the first time we present consequences of the 
fractal nature of the spectrum for the quantum time evo- 
lution of open systems. In particular, we ask how they 



are encoded in the survival probability P(t), which is the 
simplest quantity measured in laboratory experiments. 
We show that P(t) decays as 



Pit) ~ l/t a : 



a = 1 - D, 



■ 



(1) 



where Dq is the fractal (box-counting) dimension of the 
spectrum of the closed system. Moreover, we determine 
the time scale t* up to which this power-law decay can 
be observed. It scales as 



t* 



L 1 ; 7 = !/£><? 



(2) 



where L is the sample size. Beyond this time scale P{t) 
decays exponentially. Our results ([!]), (]|) are confirmed 
numerically for two types of quasi-periodic tight-binding 
models and are supported by analytical arguments. 
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FIG. 1. The survival probability P(t) of the Harper model 
(A = 2), for three different sample lengths L = 250,500,4000 
exhibits an inverse power-law P(t) ~ t~ a . A least squares fit 
yields a = 0.55 ± 0.05 in accordance with Dq ~ 0.5 and Eq. (0). 

The mathematical model we consider is the time- 
dependent Schrodinger equation 

. dip n {i) 



dt 



V n ip n {t) + ip n +i{t) + ^„-i(t), (3) 



where ipn(t) is the probability amplitude for an elec- 
tron to be at site n of a one-dimensional (ID) sample 



1 



of length L. The on-site potential V n is determined by 
quasi-periodic sequences. We assume absorbing bound- 
ary conditions at the ends of the sample |n| and initial 
excitations in the form of a 5— like packet launched at one 
of the boundaries, i.e. ip n (t = 0) = 5„ r i. Equation (||) 
has been integrated numerically using a Cayley scheme 
PH with integration time step dt = 0.1. We attached 
15 additional sites at the ends of the sample and erased 
all components of the wave packet on these sites after 
each time step dt. The decay of the norm of the wave 
packet obtained in this way was not affected by a further 
decrease of dt. 
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FIG. 2. (a) Survival probabilities P(t) of the Fibonacci model, 
for three different potential strengths Vi = 0.5, Vi = 0.75 and 
V3 = 1.5 showing inverse power-laws P(t) ~ t~ a (dashed lines). 
The sample size is L = 2000 in all cases, (b) Power-law exponents 
a of the survival probability obtained numerically as a function of 
the potential strength V for the Fibonacci model. The solid line is 
the theoretical prediction a = 1 — Dq . 

We motivate and numerically verify our results using 
first the well known Harper model, a paradigm of quasi- 
periodic ID systems with a metal-insulator transition 
[^6|,Q. It is described by a tight-binding Hamiltonian 
with an on-site potential given by 



V n = A cos(27rern) 



(4) 



This system effectively describes a particle in a two- 
dimensional periodic potential in a uniform magnetic 
field with a — a 2 eB/hc being the number of flux quanta 
in a unit cell of area a 2 . When a is an irrational number, 
the period of the effective potential V n is incommensurate 
with the lattice period. The states of the corresponding 
closed system are extended when A < 2, and the spec- 
trum consists of bands (ballistic regime). For A > 2 the 



spectrum is point-like and all states are exponentially 
localized (localized regime). The most interesting case 
corresponds to A = 2 of the metal-insulator transition. 
At this point, the spectrum is a Cantor set with fractal 
dimension Dq < 0.5 |2j| while self-similar fluctuations 
of the eigenstates exist on all scales JlGlEll . 
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FIG. 3. Resonance widths as a function of the overlapping 
elements |cfc| 2 . The data are obtained by direct diagonalization of 
the effective Hamiltonian "Heff (see Eq. (0)). The dashed lines of 
slope 1 are shown for comparison demonstrating a linear relation 
for small P In all cases the sample size is L = 1597, corresponding 
to an approximant of the golden mean a = Yg§f ■ ( a ) Harper model 
for A = 2; (b) Fibonacci model for V = 0.1; (c) Fibonacci model 
for V = 0.5; and (d) Fibonacci model for V = 1.5. 

We investigate the survival probability P(t) for the 
Harper model at the critical point A = 2. In our cal- 
culations we assume a equal to the golden mean gq — 
(\/5 + l)/2. For this case it is known that D® w 0.5 
p3| . The results for various sample lengths L are shown 
in Fig. 1. In all cases the survival probability clearly 
displays an inverse power law 



p(t) = £;i^(t)i 2 ~r 



(•5) 



The best fit to the numerical data yields a = 0.55 ± 0.05 
in accordance with Eq. (Q). 

For a further test of the validity of Eq. ([!]) we now turn 
to the Fibonacci model, where Dq can be varied. Here 
the potential V n only takes on two values ±V (V 7^ 0) 
that are arranged in a Fibonacci sequence pJ. It was 
shown that the spectrum is a Cantor set with fractal 
dimension Dq depending on V jl9]. In Fig. 2(a) we re- 
port some of our numerical results for P(t). Again we 
find a power-law decay P(t) ~ t~ a , where the exponent 
depends on the potential strength V. The exponents 
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a extracted for various V are compared with the corre- 
sponding fractal dimension Dq in Fig. 2(b) and confirm 
the validity of Eq. (TO). 



10" 

1-|C(X)| 



10 











V=0.3 






0.27 _ , ^v***/* ' 






. V=0.6_/^ 









10 



10 



10 



FIG. 4. The scattering autocorrelation function Eq. (ra) (plotted 
as 1 — C (x) I ) for some representative V values of the Fibonacci 
model compared to the theoretical expectation (dashed line). In 
these calculations the sample has length L = 10946 and is at- 
tached to one lead. The phase of the corresponding scattering 
matrix S(E) = e l *( B ) was calculated with the help of an iteration 
relation developed in Ref . Jifij . 

We now want to give a general argument for the va- 
lidity of Eq. (|TJ) . The effect of an open edge for the sys- 
tem described by Eq. (||) can be simulated by adding 
the imaginary shift i to the first diagonal element of the 
Hamiltonian matrix pi[ . Denoting our quasi-periodic 
tight-binding Hamiltonian by Hl, this approach yields 
an effective Hamiltonian 



Keff - H L 



ie, 



(0) 



where e= (1, 0, . . . , 0) T is an L— dimensional vector that 
describes at which site we impose the absorbing bound- 
ary condition. The eigenenergies of the effective Hamil- 
tonian are complex Eu—Eu — iTk/2 leading to the decay 
of the survival probability P(t). When the on-site poten- 
tial fulfills | V n | > 1 the imaginary shift can be considered 
as a small perturbation of the Hamiltonian of the closed 
system. In this case according to the perturbation the- 
ory Tfc ~ \ipi\ 2 holds, where tp* is an eigenstate of the 
closed system with energy E^. The survival probabil- 
ity is then given by P(t) ~ J2 k |cfc| 2 e _rfct , where cu are 
overlapping elements of the initial state with the eigen- 
states t/£. Choosing the initial state to be concentrated 
at site n = 1, we have |cfc| 2 = l^i | 2 ~ IV Using this and 
converting the sum into an integral we obtain: 



Pit) 



J r V{T) e - rt dT, 



(7) 



where V(T) is the resonance width distribution. In 
Ref. |H it was shown that V(T) ~ r-( 1+£> ?) for small I\ 
Inserting this expression into Eq. (j^) one finds the asymp- 
totic power-law decay stated in Eq. (m) . In order to check 



the validity of the perturbative arguments we numerically 
calculate the resonance widths and the overlapping el- 
ements Cfc for the Harper and Fibonacci models. As can 
be seen from Fig. 3 the prediction \ck\ 2 ~ of the per- 
turbation theory holds even for small V. 
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FIG. 5. Dependence of the break-time t* , on the system size L 

and the box-counting dimension Dq . (o) refer to the Harper 

model at A = 2 and various L's while (O) refer to the Fibonacci 

model with various V"s and L = 2000. The dashed line has slope 

1 and corresponds to the theoretical expectation Eq. (^). 

An immediate consequence of Eq. (|j) is that the 
scattering matrix autocorrelation function C(x) =< 
S(EyS(E+x) >e, decays in the form of a power law. In 
particular, using the relation between the survival prob- 
ability P(t) and C(x) fl3| |, we obtain 

1 - C( X ) ~ X J dtP(t)exp(-it X ) ~ X Q , X « 1. (8) 

Equation (||) is in contrast to the Lorentzian form of 
C(x) predicted by RMT for chaotic/ballistic systems 
p7[ . Comparison of C(x) for various V- values of the Fi- 
bonacci model with the theoretical prediction Eq. (||) in 
Fig. 4 shows a nice agreement and provides an additional 
check for the validity of Eq. ([!]). 

For finite samples the power-law decay of P(t) (Eq. (d)) 
holds up to a break time t* , beyond which it turns 
into an exponential decay. The rate of the latter is 
determined by the smallest resonance width T m i n and 
thus t* ~ l/r TO j n . An estimation for t* can be de- 
rived as follows. Imposing the normalization condition 
for the resonance width distribution and assuming that 
the power law V(T) ~ T- ( - 1+D o '> is valid for T > T min , 
i.e J r °° dTV(T) = L, we obtain Eq. (g). This predic- 
tion is verified numerically in Fig. 5 where we defined t* 
as the time where Pit) deviates by 5% from the power- 
law decay. We want to point out that the increase of t* 
for decreasing Dq is consistent with the enlargement of 
the interval where |cfc| 2 ~ holds and its shift towards 
smaller values (notice the change of the axes scales in 
Fig. 3(b)-(d)). 
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In summary we find that systems with a fractal spec- 
trum show a power-law decay of the survival probability 
P(t) ~ < _ ( 1_£) o ). For the finite systems of size L, this 
decay can be observed up to a time scale t* ~ L 1 / " 
beyond which an exponential decay sets in. Our predic- 
tions should be observable in propagation experiments 
with electrons or classical waves in quasi-periodic super- 
lattices or dielectric multilayers @). Moreover the con- 
nection between a, 7 and Dq provides a new possibility 
for experimental measurements of the fractal dimension 
Dq by studying current-relaxation processes. 

We thank L. Hufnagel, for useful discussions. (T.K) 
thanks U. Smilansky for initiating his interest in quan- 
tum scattering. 
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